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Abstract 

It is shown that a paradigm of classical statistical mechanics — the ther- 
malization of a Brownian particle — has a low-dimensional, deterministic 
analogue: when a heavy, slow system is coupled to fast deterministic chaos, 
the resultant forces drive the slow degrees of freedom toward a state of statis- 
tical equilibrium with the fast degrees. This illustrates how concepts useful 
in statistical mechanics may apply in situations where low-dimensional chaos 
exists. 

PACS numbers: 05.45.+b, 05.40.-hj 
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Since the study of chaotic dynamics has clarified fundamental issues in classical statistical 
mechanics it is worthwhile to consider the converse: when does intuition from statistical 
mechanics carry over to low- dimensional chaos? We all know, for instance, that a heavy 
particle immersed in a heat bath — a Brownian particle — is subject to both an average 
frictional force, and stochastic fiuctuations around this average, and that the balance be- 
tween these two thermalizes the particle. Now suppose the "Brownian" particle is coupled 
to a fast, low-dimensional, chaotic trajectory, rather than to a true heat bath. It is known 
that the particle then feels a dissipative force does the particle also (in some sense 

yet to be defined) "thermalize" with the chaotic trajectory? That is, does the fast chaos 
behave as a kind of "miniature heat reservoir" , exchanging energy with the particle in a way 
that brings the two into statistical equilibrium? In this paper, we pursue this question by 
considering the reaction forces acting on a heavy, slow system (our Brownian particle) due 
to its coupling to a light, fast trajectory. When the fast motion is chaotic, the forces on the 
particle include a conservative force, and two velocity-dependent forces, one magnetic-like, 
the other dissipative However (as in the case of coupling to a true thermal bath), there 
also exists a rapidly fiuctuating, effectively stochastic force, which has not been studied in 
detail. We describe an approach which incorporates this force, with the others, into a unified 
framework. It is shown that the inclusion of this stochastic force — related to the frictional 
force by a fiuctuation-dissipation relation — causes the slow Brownian particle and the 
fast chaotic trajectory to evolve toward statistical equilibrium. 

This result provides some justification for applying statistical arguments (involving, e.g., 
relaxation toward equipartition of energy) to physical situations of only a few degrees of 
freedom. A discussion of examples — including one-body dissipation in nuclear dynamics 
0, the Fermi mechanism of cosmic ray acceleration 0, and the diffusive transport of comets 
— where such "thermal" arguments may provide insight into the physics behind more 
explicit calculations, will be presented in Ref. 0. 

As a starting point for our discussion, we consider the framework of Ref. where the 
position R of the slow particle parametrizes the Hamiltonian h governing the fast motion: 
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h = /i(z,R), where z denotes the fast phase space coordinates. (The nature of the fast 
system will remain unspecified, but we take it to have a few, ~ 2, degrees of freedom.) 
This classical version of the Born-Oppenheimer framework has received considerable interest 
in recent years 0,§J^. We assume that, if R were held fixed, then a fast trajectory evolving 
under h would ergodically and chaotically explore its energy shell (surface of constant h) in 
the fast phase space. This sets a fast time scale, Tf, which we may take to be the Lyapunov 
time associated with the fast chaos. A slow time scale, r^, is set by the motion of the slow 
particle: it is the time required for the Hamiltonian h to change significantly. We assume 
Tf <^ Tg] thus, the fast trajectory z(t) evolves under a slowly time-dependent Hamiltonian 
h. The full Hamiltonian for the combined system of slow and fast degrees is given by 
if(R, P, z) = P^/2M + h{z, R), where P is the momentum of the slow particle, and M is 
its mass. (R, P, z) thus specifies a point in the full phase space of slow and fast variables. 
It is assumed that surfaces of constant H are bounded in the full phase space. 

Given this formulation, the force on the slow particle is F(t) = —dh/dH, evaluated 
along the trajectory z{t). From the point of view of the slow particle, this force fluctuates 
rapidly, so it is natural to separate F(t) into a slowly- changing average component, and 
rapid fluctuations F(t) around this average. In Ref. [Q, Berry and Robbins introduce an 
approximation scheme for obtaining the net average reaction force. At leading (zeroth) order 
of approximation, the ergodic adiabatic invariant p| dictates the energy of the fast system 
as a function of the slow coordinates, and this energy in turn serves as a potential for the 
slow system, giving rise to a conservative "Born-Oppenheimer" force Fq. At next order, 
the Berry-Robbins framework yields two velocity-dependent reaction forces: deterministic 
friction (Fdf) and geometric magnetism (Fgm) [|10|- Geometric magnetism is a gauge force 
related to the geometric phase; deterministic friction (see also Ref. 0) describes the irre- 
versible flow of energy from the slow to the fast variables. Thus, while at leading order the 
fast degrees of freedom create a potential well for the slow degrees, at first order the fast 
motion effectively adds a magnetic field, and drains the slow system of its energy. 

What about the effets of the rapidly fluctuating component, F(t)? If the analogy with 
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ordinary Brownian motion is correct and some sort of statistical equilibration occurs, then 
F{t) ought to play a central role in the process. We now describe a framework which 
incorporates the effects of F(t) into a description of the slow particle's evolution. 

In our framework we consider an ensemble of systems. Each member of the ensemble 
consists of a single slow particle coupled to a single fast trajectory, and represents one 
possible realization of the combined system of slow and fast variables. Representing this 
ensemble by a density in the full phase space, Liouville's equation is: 

d(j) P d(p dh d(p , r , , . ^ /.s 

where {■, ■} denotes the Poisson bracket with respect to the fast variables, z. Henceforth, we 
will ignore all information about the fast trajectory except its energy, E{t) = h[z{t),'R.{t)] 
(which evolves on the slow time scale 0). Thus, what we are really after is the evolution of 
W^(R, P, E, t), the distribution of our ensemble in the reduced space where all fast variables 
other than E have been projected out. In this reduced space, F{t) is stochastic, which in 
turn suggests that W evolves diffusively. 

The derivation of an evolution equation for W is somewhat involved, and is sketched in 
the Appendix. Here we simply state the result: 



dW P dW ^ , ,,,, e;^ 
+ D • (uW) + - A 



T,LijDj[ ^ 



(2) 



dt M dR ^ ' 2 

(Summation over repeated indices is implied.) To explain notation, we first define VL{E, R) = 

/rfz^[E-/i(z,R)]. Then S = and u = -(l/S)(9^]/9R). Next, 

^_ d _P d 
^ = dP~ MdE- 

Lij{E,TV) is an integrated correlation function defined by Eq.|20|. Finally, e ~ r/Zr^ <^ 1 is 
an ordering parameter; Eq.^ is valid to 0(e). 

Q, S, and u have simple interpretations in terms of the energy shell {E, R) [the surface 
in z-space defined by h{z, H) = E]. Q{E, R) is the volume of z-space enclosed by this shell. 
S(£',R) = / dz6{E — h) represents the statistical weight of the shell — i.e. the amount of 
fast phase space occupied by this shell — and is useful for evaluating energy shell averages: 



where {Q)e,r. denotes the average of Q{z) over the energy shell (i?,R). Finally, u(£', R) = 

{dh/dR)E,R. 

What does Eq.|^ reveal about the reaction forces on the slow particle? Below, we outline 
calculations behind the following assertions regarding the content of Eq.^: (1) it reproduces 
the average reaction forces Fq, F^/, and F^^; (2) it describes the effects of the rapidly 
fluctuating force F{t); and (3) it predicts that the Brownian particle does indeed thermalize 
with the fast trajectory. For a more detailed treatment of this problem, see Ref. [H. 

First, letting S = P^/2M + E denote the total energy of the system, note that t)S = 0. 
Thus, D is a constrained derivative: D = {d/dP)s, where the notation indicates that S, not 
E, is held fixed. This means that the evolution depicted by Eq.|^ takes place along surfaces 
of constant S in (R, P, £')-space; this is simply a statement of energy conservation. 

Next, if we explicitly separate drift terms from diffusion terms, Eq.|^ becomes 

Here, the derivatives w.r.t. P are the constrained derivatives {d/dP)£, and 

f=-u(B,R)-eA-.-^, (6) 

In Eq.^, f plays the role of a drift coefficient for the slow momentum, and thus represents the 
average force acting on the Brownian particle. A comparison with Ref. 0], Sec. 2, reveals 
that the first term of f is the leading (Born-Oppenheimer) force Fq; the second is a sum 
of the two velocity-dependent forces, F^f and F^^: if we express the matrix K in Eq.|^ as 
the sum of its symmetric and anti-symmetric components, then the former gives us F^^j, 
the latter F^.^. Eq.^ thus reproduces the average forces Fq, F^j, and Fg^ acting on the 
Brownian particle. 

The last term in Eq.^ describes the diffusion of slow momenta due to the fluctuating force 
F{t). The diffusion coefficient is the matrix L, or more precisely its symmetric component 



2^sym gy Eq.|^, however, L'^v"^ is related to K^^"^, which as mentioned is responsible for the 
dissipative force acting on the slow particle. Eq.|^ thus emerges as a fluctuation-dissipation 
relation. This relation was first noted by Berry and Robbins 0. 

Finally, do the forces acting on the slow Brownian particle cause it to thermalize with 
the fast trajectory? To answer, we must first define what we mean by "thermalization" in 
the context of the present problem (where temperature plays no role). In ordinary statistical 
mechanics, thermalization means, fundamentally, a statistical sharing of the total energy X: 
after a Brownian particle has long been in contact with a heat bath, the probability for 
finding it in some state of energy x is simply proportional to the amount of phase space 
available for the bath to have the remaining energy, X — x. (This leads to the Boltzmann 
factor P oc exp —x/ksT.) Similarly in the present context, we take the "thermalization" 
of slow and fast degrees to mean a statistical sharing of the total energy S: the slow and 
fast variables have thermalized, if the probability for finding the former in a state (R, P) 
is simply proportional to the amount of phase space available for the latter to have energy 
E = S — P'^ /2M, namely S(i?, R). (Thus, S(i?, R) plays the role of the Boltzmann constant 
here.) For our ensemble, this condition implies that an initial distribution W(R,P, E,tQ) 
evolves toward one that has the form 

W{R,P,E,t^) = Gi£)E{E,R). (8) 

[G{8) is determined by the distribution of total energies, ri{S), which remains constant.] 

We now make some formal arguments to show that Eq.^ indeed represents the ultimate 
fate of a distribution W evolving under Eq.|[ Consider an entropy S[W] = — J WlniW/T,), 
where / = J d^R J d^P J dE. Using the identity dE/dR = -{d/dE){Eu), Eq.| gives 

f^i/i;w,>o. (9) 

where T = t){W/T,). (The inequality follows from the fact that the eigenvalues of L'^y^ 
are non- negative. A proof of the latter is given in Ref. less formally, recall that the 
eigenvalues of L^y"^ are diffusion coefficients, and as such have no business being negative.) 
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Now, the distribution of total energies, fjiS), is conserved as W evolves with time. However, 
within the set of all densities W corresponding to a particular r]{S), S[W] is bounded from 
above 0. Thus as W evolves with time, the value of S[W] never exceeds a certain upper 



limit. Since dS/dt > 0, the entropy must eventually saturate, i.e. r^Oast^oo[|TT 
This in turn implies that 

W(R, P, E, t) g{S, R, t) S(E, R). (10) 

However, Eq.|lO| is a solution of Eq.^ only if g is independent of both R and t, so we 
finally conclude that W — > G{S)Ti{E,lV) asymptotically with time. Thus, the ensemble 
thermalizes, in the sense defined in the previous paragraph; this is our central result. 

This result may be restated as follows . If we start with a fast chaotic, ergodic system, 
which we then enlarge by coupling a few slow degrees of freedom to the fast ones, then the 
combined system is itself ergodic (over the surface of constant H) in the enlarged phase 
space. Thus the property of ergodicity is promoted, from the fast phase space, to the full 
phase space of slow and fast variables. 

Note also that this "thermalization" proceeds on a time scale much longer than that 
characterizing the chaotic evolution {tj). This is again similar to the case of ordinary 



Brownian motion — where such a separation of time scales is central [T^ — but stands 
in contrast to the more familiar examples of low- dimensional chaos (e.g. the N = 2 Sinai 
billiard [0), where the mixing time and the Lyapunov time are comparable. 

It is no new thing to say that a chaotic, ergodic trajectory offers a low-dimensional 
(A^ ~ 2) analogue for a truly thermal (A^ ^ 1) system. The novelty of the present work 
is that it extends this analogy to encompass the important paradigm of Brownian motion, 
where the thermal system or chaotic trajectory is coupled to a few degrees of freedom 
characterized by a much longer time scale. Then, in either case, the forces acting on the 
slow system drive it toward a state of genuine statistical equilibrium with its environment. 

Finally, it would be interesting to study the quantal version of this problem. Srednicki 
1^ has recently argued that concepts from quantum chaos may provide a solid foundation 



for quantum statistical mechanics. The focus in Ref. [0 is on genuinely thermal systems 
(A^ ^ 1), and does not deal specifically with the case when a few degrees of freedom are 
slower than the rest. Nevertheless, Srednicki's application of Berry's conjecture |[T^ to the 
quantal evolution of a classically chaotic system might serve as a guide to a quantal analysis 
of the purely classical problem studied here. (To the best of my knowledge, no one has looked 
explicitly at the application of Berry's conjecture to a system which classically exhibits two 
widely separated time scales.) 

It is a pleasure to acknowledge that conversations with Greg Flynn, Allan Kaufman, 
Robert Littlejohn, Jim Morehead, and Wladek Swi§tecki were very useful in obtaining the 
results presented in this papers. This work has been supported by the Department of Energy 
under Grants No. DE-AC03-76SF00098 and DE-FG06-90ER40561, and by the National 
Science Foundation under Grant No. NSF-PYI-84-51276. 

Appendix 

Here we sketch the derivation of Eq.|^ from Eq.^ using what is essentially (though not 
explicitly) a multiple-time-scale analysis, and is similar to that of Ref. [Q]. 

To begin, we use our adiabaticity parameter e -C 1 to formally incorporate into Eq.^ the 
assumption that (R, P) is "slow and heavy", whereas z is "fast and light": 

90 P 90 dh <90 . 

^a^ + ^M-9R-^M-9P + ^'^'^> = °- 

With this modification, changes in (R, P) take place over times of order unity, whereas 
changes in z occur over times of order e. 

Next, we insert the Ansatz = 0o -|- e0i + e^02 into Eq.|ll] and order by powers of 

e: 

{/i,0o} = O (12) 

r, . 1 f d P d dh 9 \ 
l^'M=(a^ + M-M-9R-9pj'^-^' 

r > 1. Since h commutes only with functions of itself under the Poisson bracket (by assump- 
tion of ergodicity), the solution to Eq.|I^ has the form 0o(R, P,z,t) = A(Tl,P,h{z,'R),tj. 
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To solve for the dependence of A on its arguments, we must examine Eq.|T3|, with r = 1: 

Taking a phase space average of both sides over some energy shell {E, R) in the fast phase 
space, we get 

dA P OA ^ ^ 

where the third argument of A is now E, the energy of the shell over which the average is 
taken. With the identity D ■ (Su) = (P/M) ■ 9S/(9R, we rewrite Eq.|r3| as 

To solve for we first use Eq.|15| to rewrite Eq.p^: 

adopting the notation of Ref. 0]. [The left side of this equation is evaluated at (R, P, z, t); 
the value of the third argument of A on the right side is E = /i(z,R).] This has the 
form {h, /} = g; the general solution 0] consists of both a homogeneous term, (pin = 
5(R, P, h, t), and an inhomogeneous term, 

/■o dh 

0u(R, P, z, t) = ds — (z„ R) ■ DA, (18) 

where Zs(z,R) is the point in phase space reached by evolving a trajectory from z, for a 
time s, under the Hamiltonian /i(z,R). Note that {4>ii)e,'r = for any (i?,R). 

We solve for B much as we did for A: writing Eq.|l^, with r = 2, we average each side 
over an energy shell (-E',R). After manipulation, this gives 



where 



fO / dh f dh\ \ 
^^. = 2/ ds{-- --) . (20) 



— oo 



\ ' \ J / SI i?.R 



The first factor inside angular brackets is evaluated at z, the second at Zs(z, R); the average 
is over all points z on the energy shell {E, R). It is assumed that the integral converges. 
Finally, iy(R, P, E, t) is given by a projection of from (R, P, z) to (R, P, E): 



W 



J dz6{E-h)<f) = J:{E,R) (0)^,K. 



(21) 



Since {4>)e,'r = A + eB (to order e), we combine Eqs.|T6|, and ^ to obtain the desired 
result, Eq.||. 
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